Random Antiferromagnetic Spin-1/2 Chains with Competing Interactions 
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We study disordered antiferromagnetic spin-1/2 chains with nearest- and further-neighbor in- 
teractions using the real-space renormalization-group method. We find that the system supports 
two different phases, depending on the ratio of the strength between nearest-neighbor and further- 
neighbor interactions as well the bond randomness strength. For weak further neighbor coupling 
the system is in the familiar random singlet phase, while stronger further neighbor coupling drives 
the system to a large spin phase similar to that found in the study of random antiferromagnetic- 
ferromagnetic spin chains. The appearance of the large spin phase in the absence of ferromagnetic 
coupling is due to the frustration introduced by further neighboring couplings, and is unique to the 
disordered chains. 
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I. INTRODUCTION 

One-dimensional (ID) quantum spin systems have 
been of interest to physicists for many years. This is 
not only because these systems have been good testing 
grounds for various theoretical techniques and approxi- 
mations but also because they exhibit a wealth of fas- 
cinating low-energy physics. Among various intriguing 
phenomena of these systems, the interplay between quan- 
tum fluctuation and disorder has attracted considerable 
recent attention. The most thoroughly studied model in 
this context is the random antiferromagnetic (AF) spin- 
1/2 chain with nearest neighbor interaction. It has been 
shown/* using the celebrated real space renormalization 
group (RSRG) method^ that the low-energy physics 
of the model is controlled by the random singlet (RS) 
fixed point of the RSRG and is universal. Among the 
universal properties of the random singlet phase are the 
uniform spin susceptibility: \ ~ l/Tlog 2 T, and the dis- 
order averaged spin-spin correlation function: (Si ■ Sj) ~ 
(-iy- j /(i-j) 2 . The RSRG method (with proper exten- 
sions) has also been applied with considerable success to 
a number of other disordered spin chain models (all with 
nearest neighbor interaction only) jiS&LSuLiSiiiii^ii 3 . as 
well as two-leg spin ladders pliiiiiiii 

In the present work we study random AF spin-1/2 
chains with nearest and further neighbor couplings, using 
the RSRG method. Our motivation comes from the fol- 
lowing considerations. First of all, as mentioned above, 
existing theoretical studies have been focusing on mod- 
els with nearest neighbor couplings only; the RG flow 
equations of the couplings are relatively simple in this 
case which allows, for example, exact analytical solution 
of the fixed point in the case of random AF spin-1/2 
chain. 1 In real physical systems, on the other hand, fur- 
ther neighbor couplings are always present, and in cer- 
tain cases they can even be quite strong. There are a few 
promising experimental realizations of materials that ex- 
hibit non trivial next-nearest neighbor interactions. One 
of the examples of real physical systems that may meet 



the criteria is CuGeQ3i 17 i 18 ' 19 i 20 ' 21 Studies on this system 
have revealed that the angle of Cu-O-Cu bond is close to 
90°. This will induce a competition between antiferro- 
magnetic superexchange between the Cu ions mediated 
by the oxygen ion and ferromagnetic direct exchange be- 
tween the Cu ions. As a result the nearest-neighbor su- 
perexchange interaction is weakened and hence it is ex- 
pected that the next-nearest-neighbor interactions which 
arise from the Cu-O-O-Cu path cannot be neglected. 
The strength of the second-neighbor bonds can also be 
controlled by applying pressure to such systems. Ma- 
suda and coworker— studied the effect of pressure on 
highly Mg-doped CuGeOa and found that the frustra- 
tion is enhanced as the pressure is increased. Another 
example of material that exhibits non trivial second- 
neighbor interaction is CugGe6 Ois — xH^O studied by 
Hase and coworker. 23 Motivated by these experimental 
realizations, we would thus like to study the effects of 
next-nearest-neighbor interactions, and in particular, the 
stability of the RS fixed point against their presence. 

Secondly, nearest neighbor models have no frustration 
in them. Further neighbor interactions, on the other 
hand, can introduce frustration, and this is known to lead 
to new physics and phases in the case of pure chains. For 
example, it is known in the case of spin-1/2 chain with 
nearest and next nearest neighbor couplings (Ji and J%), 
there are two different phases depending on the ratio be- 
tween the twoi 24 i 25 i 26 i 27 i 28 For zero or small J2/J1, the 
system is in a gapless (critical) phase with power-law 
spin-spin correlation, while for larger J2/J1 the system 
spontaneously dimerizes and opens a gap in the exci- 
tation spectrum, and the spin-spin correlation becomes 
short-ranged. In the special case of J2/J1 = 1/2, which is 
the so-called Majumdar- Ghosh model, the ground state 
of the system is known exactly; they are collections of 
neighboring spins forming singlet pairs over either even 
or odd nearest neighbor bonds ^^SL^ It is thus of in- 
terest to study how frustration affects the physics of dis- 
ordered chains, and whether new phases can be stabilized 
by it. 
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Our results can be summarized as follows. We find 
that there are two phases in the model we are consid- 
ering, controlled by the ratio of the strength of nearest- 
neighbor and next-nearest neighbor interactions and the 
strength of bond randomness. The Random Singlet (RS) 
phase is found to be stable against weak further neighbor 
couplings; in this case the strength of further neighbor 
couplings (as measured by strength of nearest neighbor 
couplings) flow to zero as energy scale decreases, thus 
the low-temperature properties of the system is still con- 
trolled by the Random Singlet (RS) fixed point. For 
strong enough further neighbor couplings, on the other 
hand, the RS phase becomes unstable and the system is 
driven into another phase which is controlled by large 
effective spins at low energies. We find that in this 
phase the system is still dominated by effective nearest- 
neighbor interactions at low-energy; however the effec- 
tive couplings can be either antiferromagnetic or ferro- 
magnetic, with random distributions. We conclude that 
this phase is the same as that found in random antifer- 
romagnetic (AF)-ferromagnetic (F) spin chain systems 
with nearest-neighbor interactions only studied by West- 
erberg et al.A The physical origin of the appearance of 
effective ferromagnetic couplings is the frustration intro- 
duced by further neighbor couplings. 

The remainder of the paper is organized as follows. In 
Sec. II we introduce the model we study and discuss 
the application of RSRG to this model. Results of our 
numerical studies on the model are presented in Sec. III. 
In Sec. IV we summarize our findings and make contact 
with previous works that are related to our studies. 

II. THE MODEL 

We consider the antiferromagnetic (AF) spin-1/2 chain 
described by the following Hamiltonian : 

JV-l N-2 
i=l i=l 

where N is the number of spins on the chain, Si is a 
spin-1/2 operator at the i th site and the positive cou- 
plings Ji and Ki are distributed randomly according to 
some probability distributions which will be described 
in more details in the next section. The Hamiltonian 
written down in Eq. consists of two terms, where 
the first term describes nearest-neighbor interactions be- 
tween the spins and the second term describes next- 
nearest-neighbor (n.n.n.) interactions. The schematic 
diagram of the system described by the Hamiltonian 
is depicted in Fig. [TJa). We will mostly focus on chains 
with n.n. and n.n.n couplings in this paper; but some 
results of chains with couplings beyond n.n.n. will also 
be presented. 

We use the real-space renormalization-group method 
to study the Hamiltonian (JIJ. The application of this 
method to AF spin-1/2 chain with n.n. couplings only 



is well known. The basic idea is to isolate the strongest 
bond in the system, decimate it, and calculate the effec- 
tive interactions generated between what were the third- 
nearest neighbors. The key simplifying features in this 
case are that the generated interactions are always an- 
tiferromagnetic, and they connect only nearest-neighbor 
spins (after the two spins coupled by the strongest bond 
are removed). 

Appropriate extensions of the original RG scheme need 
to be included in order to study the present model with 
further neighbor couplings properly. First we notice that 
the coordination number, i.e. the number of spins cou- 
pled to a given spin, grows as the energy scale is low- 
ered so we need to keep track of the structure of the sys- 
tem. This is in contrast to the AF spin-1/2 chain with 
n.n. couplings where the coordination number is always 
2. Second, as we will see later, effective ferromagnetic 
couplings may be generated at certain stage as RSRG 
is carried out in the presence of antiferromagnetic n.n.n. 
couplings. The formation of ferromagnetic couplings al- 
lows the possibility of generating effective spins with sizes 
larger than 1/2, so we need to extend the RG rules to in- 
corporate arbitrary spin sizes and coupling signs. Let us 
discuss these in more detail. Consider spin 3 and 4 in 
Fig. |TJa) which are coupled by the strongest bond, and 
other spins in the system that couple to at least one of 
them. Due to the presence of n.n.n. couplings, we have 
a six-spin problem instead of a four-spin problem for a 
given pair of spins coupled by the strongest bond. The 
Hamiltonian for the six-spin problem is given by : 

H = H + Hj, (2) 

where 

Ho = • S4, 

Hi = J23S2 ■ S3 + J45S4 ■ S5 + J13S1 ■ S3 

+ J35S3 • S5 + J24S2 ■ S4 + J4QS4 ■ S4, (3) 

where Jij is the antiferromagnetic coupling between Si 
and Sj . We have shown in our previous work on the spin 
ladders^ that to the second order perturbation calcula- 
tion, Hi only generates pairwise interactions among the 
spins and hence it is only necessary to include a pair of 
spins coupled to the two spins connected by the strongest 
bond when we consider the effective interaction between 
them, i.e. we just have to consider four-spin clusters for 
a given segment which contains the strongest bond. Let 
us consider the most complicated four-spin cluster where 
a given spin is coupled to three other spins as depicted 
in Fig. |3 The renormalized coupling between two spins 
in the cluster, say spins 2 and 5, is given by : 

- 1 

J25 — J25 + 777 — {J23J45 + J24J35 — J23J35 — J24J45) 
ZJ34 

= J25 + 7n~( J 23 - J24XJ45 - ^35) (4) 

ZJ34 

where Jy is the renormalized coupling between Si and 
Sj, and is the original bond between Si and Sj. 
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Examining Eq. we can see that some of the contribu- 
tions to the renormalizcd coupling from 2nd order pro- 
cesses are ferromagnetic. The overall sign of the total 
interaction between the second and fifth spin will be de- 
termined by the relative strength between the antifer- 
romagnetic nearest-neighbor and next-nearest-neighbor 
bonds. In general if the n.n.n. couplings are very weak 
compared to the n.n. couplings then the ferromagnetic 
interactions will not appear. This is quite different from 
what we found in the study of the ladder where effec- 
tive ferromagnetic interactions appear as soon as the RG 
is applied to the system. Due to the possibility of the 
appearance of ferromagnetic couplings at some step of 
RG, it is necessary to generalize the RG procedure to 
include arbitrary spin sizes and coupling signs. The dis- 
cussion on how this is done has been spelled out in great 
detail in our earlier work on the spin ladder^ We carry 
out the numerical calculation using the rules described in 
previous paragraphs and present the results in the next 
section. 



III. NUMERICAL RESULTS 

We present numerical results for spin chains with 
nearest-neighbor (n.n.) and next-nearest-neighbor 
(n.n.n.) interactions with total number of spins up to 
60000. We search for the bond with the largest gap, A , 
which is defined as the gap between the ground state 
and the first excited state, decimate it, and calculate the 
effective interactions among the remaining spins. The 
procedure is repeated until the number of spins left is 
about 1% of the original number of spins in the system. 
We use 100 samples and take the disorder average over 
all these samples in all our calculations. The nearest- 
neighbor bonds are chosen to be distributed randomly 
according to the power-law probability distribution: 

P n . n .{Ji) = {l-a)J^\{) < Ji<l, (5) 

where the power-law exponent a < 1 parametrizes the 
randomness strength; the larger a, the stronger the ran- 
domness. The reason for choosing a power-law form is 
because for the random spin-1/2 chain, the fixed point 
distribution is known to be in the power-law form. So by 
choosing initial distributions in the power-law form, we 
expect to start closer to the fixed point and hence reduce 
the necessity to use larger system size. 

We consider two different ways of generating the n.n.n. 
bonds. First we consider n.n.n. bonds which are com- 
pletely correlated with the n.n. bonds, where next- 
nearest-neighbor bond Ki is determined from the n.n. 
bonds through the following relation: 

"0 

where A is a parameter introduced to control the strength 
of next-nearest-neighbor interactions and Qq is the cutoff 



of the initial nearest-neighbor bonds distribution, which 
is 1. In the limit A — > 0, the AF spin-1/2 chain with 
nearest-neighbor interactions only is recovered. Eq. © 
comes from the following consideration. The interactions 
between two spins come from the overlap integral of the 
electron wave functions which are bound to the atoms 
sitting on the lattice sites. In general, the wave function 
decays exponentially at large distances, and so does the 
overlap integral. Let us consider three electrons sitting on 
different lattice sites labeled 1,2, and 3. For two electrons 
separated by a distance R, the typical interaction would 
have a form of J ~ e~ R l a , where a is a length scale of 
order the size of the wave function. Based on this picture, 
the interaction between the first and third spins, which is 
basically the overlap integral between the first and third 
spins, can be written as J ~ e -{R3-Ri)/a^ wnere ^ 3 anc j 
Ri are measured with respect to some reference point. 
This relation can be rewritten as : 

J „ e '(R 3 -R2)/a e -(R 2 - Rl )/a ^ j 2 j u (7) 

where Ji is the overlap integral between Si and Si+i. 
Hence, it is reasonable to model the correlation as the 
product of two nearest-neighbor bonds as shown in Eq. 
(0. We focus mostly on this type of further neighbor 
coupling, and unless stated otherwise, the results pre- 
sented below are for this type of further neighbor cou- 
pling. For comparison, we have also studied cases in 
which the n.n.n. couplings are uncorrelated case with 
the n.n. couplings, i.e. the n.n.n. bonds are distributed 
randomly in the system, independent of of the distribu- 
tion of the n.n. bonds. We choose the distribution to 
be in a power-law form with the same exponent, but a 
different cutoff A : 

P n .n.n.(Ki) = ]^K~ a ,0 <K l <K. (8) 

Again A parametrizes the strength of n.n.n. couplings. 
As we will see later, while the topology of the phase 
diagrams are the same for these two cases, there is 
huge quantitative differences in the position of the phase 
boundary. 

As we carry out RSRG numerically, we monitor the ap- 
pearance and proliferation of large effective spins in the 
system. We plot the sample-averaged fraction of spins 
larger than 1/2 as a function of energy scale, A , in Fig. 
[3] The left panel of Fig. [3] shows how the formation 
of large effective spins evolves as the energy scale, Ao, 
is lowered by fixing a = and varying n.n.n. bond 
strength controlled by A, while the right one by fixing 
A = 0.55 and varying a. Let us analyze the left panel of 
Fig. |21 It is very clear that, for fixed a, different anti- 
ferromagnetic n.n.n. bond strength will lead to different 
scenarios in the low energy limit. For weak enough A (in 
the regime where A < 0.5) we do not find spin sizes other 
than 1/2; not only we never find any spin larger than 1/2 
but also we never find any ferromagnetic bonds in this 
regime. The situation drastically changes when we tune 
the strength of antiferromagnetic n.n.n. bonds up to 0.55 
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where we can see clearly that large effective spins dom- 
inate in the low energy limit and drive the system into 
a new phase. This can be understood in the following 
way. For weak enough n.n.n. bonds, these interactions 
are always suppressed by the presence of n.n. bonds. We 
have explained in Eq. (J3J that the ferromagnetic bond 
will appear if the n.n.n. bonds are strong enough to over- 
come the n.n. bonds. Apparently for A < 0.5, the n.n.n. 
bonds are too weak to compete with n.n. bonds so we 
never see the emergence of ferromagnetic interactions in 
the system. On the other hand, for A > 0.55, the antifer- 
romagnctic n.n.n. bonds are strong enough to overcome 
the n.n. bonds and allow the appearance of ferromag- 
netic bonds which in turn will drive the system into a 
new phase controlled by large effective spins. 

The right panel of Fig. [3] shows another study of how 
large effective spins appear in the system by varying the 
disorder strength a for fixed A = 0.55. We find that 
the formation of large effective spins is suppressed as the 
bond disorder gets stronger. This also has a simple ex- 
planation. With increasing bond disorder strength, the 
probability of finding weak n.n. bonds is getting bigger. 
This will give us even weaker n.n.n. bonds because of the 
correlation between a next-nearest-neighbor bond with 
two nearest-neighbor bonds as given by Eq. [S] These 
weak n.n.n. bonds can not compete with the n.n. bonds 
which in turn will suppress the formation of ferromag- 
netic bonds in the system. Based on this view, we can 
understand why large effective spins are more difficult to 
form in the regime where the bond disorder is strong. So 
for strong enough bond disorder, no ferromagnetic bonds 
will appear due to the fact that n.n.n. bonds could not 
compete with n.n. bonds and the system will remain in 
the Random Singlet (RS) phase. 

The appearance of a new phase can also be deduced 
from plotting the sample-averaged \F as a function of 
temperature where the temperature is associated with 
the energy scale, Aq. We plot this in Fig. 0]where in the 
left panel a is fixed and A is varied, whereas in the right 
panel A is fixed and a is varied. \F in the RS phase 
is well known to be given by \/log 2 T. For fixed a = 0, 
we can see increasing deviations from l/log 2 T with in- 
creasing strength for A which gives us a clear indication 
that the system is driven away from the RS phase; for 
A > 0.5 instead of falling as l/log 2 T, %T appears to ap- 
proach a constant in the low-T limit. The explanation 
for this behavior is similar to the discussion in the pre- 
vious paragraph. Strong enough n.n.n. bonds will allow 
the appearance of ferromagnetic bonds which in turn will 
form large effective spins in the low energy limit. These 
strongly correlated effective spins govern the susceptibil- 
ity of the system at low temperature. The susceptibility 
in this phase has different origin from the susceptibility 
for the RS phase where the contribution comes from the 
undecimated half spins. The same situation is encoun- 
tered when A is fixed and a is varied as shown on the 
right panel of Fig. 01 The deviations are more signifi- 
cant for small a. This is consistent with our discussion 



on the previous paragraph that for strong enough bond 
disorder, the system will remain in the RS phase because 
the overall strength of n.n.n. bonds is much weaker than 
that of n.n. bonds. This is indeed what we see in our 
numerical results that \F for bigger a(> 0.6) is closer to 
the value for RS phase l/log 2 T. 

We have established that there two phases in the 
system. The transition from one phase to another is 
controlled by the strength of bond disorder a and the 
strength of n.n.n. bonds A. For a = and A < 0.5 the 
system remains in the RS phase while for A > 0.55 the 
system is driven into the new phase. We have already 
seen that the new phase is controlled by large effective 
spins in the low energy limit. Is there any other pa- 
rameter we can use to study the nature of of the new 
phase? We address this question by studying the ra- 
tio of n.n. bond strength to n.n.n. bond strength in 
the two phases as shown in Fig. |SJ It is found that 
on either side of the phase boundary, nearest-neighbor 
bonds always dominate further-neighbor bonds^ Now 
we have a more complete picture of the new phase found 
in the system. The new phase is controlled by large 
effective spins in the low energy limit and the domi- 
nant interactions come from the nearest-neighbor bonds 
only. These nearest-neighbor interactions consist of both 
antiferromagnetic and ferromagnetic bonds. These re- 
sults suggest that in the low energy limit, spin chains 
with antiferromagnetic n.n. and sufficiently strong n.n.n. 
interactions behave just like random antiferromagnetic- 
ferromagnetic spin chains, including a Curie susceptibil- 
ity discussed earlier. This brings us to the conclusion 
that the new phase found in the system we are studying 
is the same as the large spin phase found in the ran- 
dom antiferromagnetic-ferromagnetic spin chains in the 
low energy limit. The numerically determined phase dia- 
gram for spin chains with random antiferromagnetic n.n. 
and n.n.n. bonds is shown in Fig. 

The left panel of this figure shows the phase diagram 
for the correlated next-nearest-neighbor bonds as given 
by Eq. 10 whereas the right panel for uncorrelated next- 
nearest- neighbor bonds. In both cases we find that the 
system supports only two phases, which are the Random 
Singlet phase and the Large Spin phase. There are some 
differences of the phase boundaries on these two cases. 
First, the trend on how the phase boundaries change as 
we vary a and A is different for the correlated and uncor- 
related next-nearest-neighbor bonds. For the correlated 
case, A stays constant as we increase a from to 0.6 and 
tends to increase for a larger than 0.6. For the uncorre- 
lated one, A decreases with increasing a. Secondly, the 
magnitude of critical A for the uncorrelated n.n.n. cou- 
plings is much smaller than that for the correlated case, 
by as much as 10 orders of magnitude for a close to 1! 

We believe the these differences can be understood as 
the following. For the uncorrelated case we assign a prob- 
ability distribution function for the n.n.n. bonds whose 
cutoff is determined by A, and the bonds are generated 
independent of the configuration of the n.n. bonds. Al- 
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though in general the strength of the n.n.n. bonds is 
much weaker than that of n.n. bonds when A is small, 
due to the uncorrelated nature of the way they are gen- 
erated, there is a small probability that the next-nearest- 
neighbor coupling is actually stronger than the nearest- 
neighbor ones in some region of the system. As we have 
explained earlier in the text, the overall sign of the total 
interaction generated by RG between two spins depends 
heavily on the relative strength of the antiferromagnetic 
n.n. and n.n.n. bonds; thus such rare events can lead 
to the generation of ferromagnetic bonds, which in turn 
may proliferate as energy scale goes down. In the corre- 
lated case, on the other hand, such rare events are greatly 
suppressed by the correlation between n.n. and n.n.n. 
bonds. We also know that a parametrizes the width of 
the distribution; for a given A in the uncorrelated case, 
the bigger a is, the wider the distributions are for both 
the n.n. and n.n.n. bond distributions, thus the larger 
the probability of the rare events discussed above are, 
and the more likely ferromagnetic couplings get gener- 
ated. On the other hand this effect is again suppressed 
for the case of correlated n.n.n. bonds, due to the way we 
parametrizes their strength; the larger a is, the smaller 
the overall strength of the n.n.n. bonds due to the way 
the are generated. 

As discussed earlier, the appearance of effective ferro- 
magnetic couplings is a consequence of competition be- 
tween nearest and further neighbor couplings, or frus- 
tration. We have also studied spin chains with further- 
neighbor interactions which do not introduce frustra- 
tion to the system. This can be done by introducing 
ferromagnetic next-nearest-neighbor bonds or antiferro- 
magnetic third-nearest-neighbor bonds. The ferromag- 
netic next-nearest-neighbor bonds and antiferromagnetic 
third-nearest-neighbor bonds are generated in the same 
way as discussed at the beginning of this section, i.e. the 
bonds are generated through Eq. iJBJ. We present our 
results for this particular system in Fig. [7] 

The upper panels of Fig. [7|show the sample-averaged 
plot of the strength of nearest-neighbor interactions com- 
pared to the strength of further-neighbor interactions and 
the fraction of spins with sizes larger than 1/2 as a func- 
tion of energy scale, Ao, for the system with ferromag- 
netic next-nearest-neighbor bonds. We choose to fix a = 
and to vary A to see how the ratio changes as the en- 
ergy scale is lowered. We find that the nearest-neighbor 
interactions always dominate over further-neighbor inter- 
actions at all energy scale. The evolution of the spin sizes 
as the energy scale is lowered is also studied here. The 
result shows that no spin having sizes larger than 1/2 
is found in the system. Based on these results we con- 
clude that the presence of ferromagnetic next-nearest- 
neighbor bonds does not drive the system into a new 
phase. The couplings are dominated by antiferromag- 
netic bonds which suppress the formation of effective 
spins larger than 1/2 at low energy. In the low energy 
limit the system stays in the RS phase. The lower panels 
of Fig. [7| show the plot of the ratio between the strength 



of nearest and further neighbor bonds and the fraction of 
spins with sizes larger than 1/2 as a function of energy 
scale, Ao, for the system with antiferromagnetic third- 
nearest- neighbor bonds. We also fixed a = and vary A 
for this case. The results are the same for those with fer- 
romagnetic next-nearest-neighbor bonds. These results 
give us a strong indication that the system stays in the 
RS phase. We can thus conclude that non-frustrating 
further neighbor bonds act as irrelevant perturbations in 
the low energy limit, and hence the system stays in the 
RS phase. 



IV. SUMMARY AND DISCUSSION 

In this paper we have used the real space renormaliza- 
tion group method to study random antiferromagnetic 
spin-1/2 chains, with both nearest- and further neigh- 
bor interactions. We find that the system supports two 
phases, the random singlet phase and the large spin 
phase. The latter is only stabilized by sufficiently strong 
further neighbor couplings that compete with the near- 
est neighbor couplings, so that there is frustration in the 
system. 

The real space renormalization group procedure is 
quantitatively accurate only when the initial distribu- 
tions of the couplings are broad. We believe, however, 
our conclusion remains valid even if the initial distri- 
bution of couplings is not broad. In the case of near- 
est neighbor coupling only, Doty and Fisher- 4 showed 
that weak bond randomness is a relevant perturbation 
that immediately destabilizes the Luttinger liquid fixed 
point that describes the gapless phase of the pure chain, 
and bond randomness grows as energy scale goes down, 
eventually brings the system to the random singlet fixed 
point. Their arguments remain valid even in the presence 
of further neighbor couplings, as long as they are not 
strong enough to destabilize the gapless phase in the ab- 
sence of bond randomness. On the other hand when they 
are strong enough to put the pure system in the gapped 
phase with spontaneous dimerization, one of us^ showed 
that the dimerized phase is also unstable against weak 
randomness, as randomness nucleates solitons and de- 
stroys spontaneous dimerization; the low-energy degrees 
of freedom are the half-spins carried by the solitons, with 
random interaction with broad distribution (due to the 
fluctuation of inter soliton distance etc). Depending on 
whether the coupling between these spins are purely AF 
or both F and AF, the systems can be in either one of 
the two phases we find here. We thus conclude these are 
the only two phases the system supports in the presence 
of any amount of bond randomness^ 

The frustration induced ferromagnetic coupling and 
the resultant large spin formation has been discussed in a 
different context before^ In that work Yang and Bhatt 
studied spin-1 chains with random AF nearest neighbor 
bonds, with both quadratic and bi-quadratic couplings 
on each bond. It was shown that even though overall 
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each individual bond is AF, as long as in some of the 
bonds the quadratic and bi-quadratic couplings have op- 
posite tendency (i.e., one AF and the other F), effective 
ferromagnetic couplings may be generated at low-energy, 
and the large spin phase stabilized. In this case the bonds 
are frustrated due to the competition between quadratic 
and bi-quadratic couplings on the same bond. Thus the 
phenomena of frustration induced large spin formation, 
although never seen in pure systems, may actually be 
rather generic in disordered systems. 
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FIG. 1: (a) Schematic diagram for the AF spin-1/2 chain given by the Hamiltonian Q. In addition to the nearest-neighbor 
couplings between the spins, we also include the next-nearest-neighbor couplings represented by the dashed lines. Here the 
strongest bond is represented by the thick bold line, (b) The renormalization scheme after the strongest bond is decimated. 
The thick dashed lines are the renormalized couplings. 
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FIG. 2: The most complicated structure of a four-spin cluster where a given spin is coupled to the other three spins. 
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FIG. 3: The sample-averaged fraction of spins larger than 1/2 as a function of energy scale, Ao. The error bars are about 
the size of the data points shown in the figure. The left panel shows how the fraction of spins larger than 1/2 for a = 0.0 
changes as A is varied and the right panel for A = 0.55 as a is varied. Both are calculated for N=60,000. Strong enough 
next-nearest-neighbor interactions will drive the system into a new phase controlled by large effective spins. All calculations 
are done with correlated next-nearest-neighbor bonds given in Eq. ©. 
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FIG. 4: The sample-averaged yT as a function of parameters of the model, a and A. The error bars are about the size of 
the data points. For strong enough correlated next-nearest-neighbor interactions, given in Eq. ©, the susceptibilities behave 
differently from l/Tlog 2 T. The contribution to the susceptibilities come from large effective spins formed at low temperature. 
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FIG. 5: The sample-averaged ratio of the strength of the nearest-neighbor bonds to the strength of the bonds that are beyond 
nearest neighbor as a function of energy scale. It is clear from the plot that in either side of the phase, the interactions are 
dominated by nearest-neighbor bonds only. We use the correlated next-nearest-neighbor interactions defined in Eq. JjJ. 
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FIG. 6: (a) The numerically determined phase diagram for spin chains with competing interactions between nearest-neighbor 
and next-nearest-neighbor interactions. The n.n.n. interactions are correlated with the n.n. interactions (see text), (b) The 
numerically determined phase diagram for spin chains with uncorrelated n.n.n. interactions. In both cases a denotes the 
strength of the bond randomness and A represents the strength of the next-nearest-neighbor interactions. The crosses in both 
figures are obtained from numerical calculations. The dashed lines are drawn by connecting the data points to see the phase 
boundary more clearly. 
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FIG. 7: The sample-averaged ratio of the strength of the nearest-neighbor bonds to the strength of the bonds that are beyond 
nearest neighbor and the fraction of spins larger than 1/2 as a function of energy scale as a function of energy scale for 
the model with no frustration introduced into the system. Two types of interactions which do not generate frustration, e.g. 
ferromagnetic second neighbors and antiferromagnetic third neighbors, are introduced into the system. The upper two panels 
show the calculation for a model in which ferromagnetic second-neighbor interactions are introduced into the system while the 
lower two panels for antiferromagnetic third-neighbor interactions. All graphs are calculated for a — but with varying A. It 
is clear from the plot that the interactions are dominated by nearest-neighbor bonds only regardless the value of A and there 
is no formation of effective spins whose sizes are larger than 1/2. 



